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RADAR  RESONANCE  REFLECTION  FROM  SETS  OF 
PLANE  DIELECTRIC  LAYERS 


P.  D.  Jackins  and  G.  C.  Gaunaurd 
Naval  Surface  Weapons  Center 
White  Oak,  Silver  Spring,  MD  20910 


ARSTRACT 


The  prediction  of  the  Resonance  Scattering  Theory  ( R ST )  for  the 
reflection  coefficient  from  a  set  of  two  contiguous  plane  dielectric 
layers  separating  two  semi-infinite  dissimilar  non-conducting  media, 
is  constructed  and  compared  to  the  exact  classical  model  solution. 
The  comparison  serves  to:  a)  show  the  accuracy  and  simplicity  of  the 
RST-predi cti on ,  and  b)  to  underline  the  usefulness  of  the  RST  to 
produce  simple  physical  interpretations  of  generally  complex 
phenomena.  The  analysis  provides  a  systematic  method  for  detecting 
the  presence  of  a  dielectric  1 ayer  under  another  one  covering  it 
(possibly  the  situation  caused  by  an  oil  spill  in  ice-covered  Arctic 
regions),  by  certain  modulation  effects  present  in  the  "response 
surface"  of  the  returned  echoes.  This  method  also  identifies  the 
material  composition  of  the  lower  or  hidden  layer  in  the  bilaminar 
configuration.  The  process  disentangles  which  resonance  feature 
present  in  the  radar  reflection  coefficient  is  caused  by  which  of  the 
two  interacting  layers.  This,  therefore,  solves  the  inverse 
scattering  problem  for  the  composi ti on  and  thickness  not  only  of  the 
top  visible  (ice)  layer,  but  also  of  the  substance  (oil)  hidden  under 
that  upper  layer. 


INTRODUCTION 


For  many  years  wave  reflection  and  transmission  through 
penetrable  layered  media  has  been  a  well  plowed  field  of  research 
where  solutions  usually  have  been  given  in  terms  of  reflection  and 
transmission  coefficients.  These  coefficients  describe  the  energy 
fraction  returned  by,  or  transmitted  through,  the  (possibly  layered) 
scatterer  (1-7).  While  such  solutions  are  exact,  they  are  so 
cumbersome  for  systems  of  two  or  more  plane  Interacting  layers  that 
they  hardly  ever  permit  a  quick  understanding  of  the  physics  behind 
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the  results.  A  novel  analytical  approach  denoted  the  Resonance 
Scattering  Theory  ( R ST )  (8)  has  been  recently  developed  at  NSWC. 
This  method  has  been  applied  to  many  probelms  dealinp  with  the 
scattering  of  (acoustic,  el ectromagnetic  and  elastic)  waves  by 
penetrable  objects,  mostly  of  simple  shapes  (8).  While  the  results 
generated  by  means  of  the  RST  are  mostly  approximations,  their 
simplicity  and  accuracy  permits  a  greater  Insight  and  deeper 
understanding  of  the  physical  mechanisms  involved  in  their  generation 
than  was  ever  before  possible. 

In  qeneral,  under  the  dissecting  knife  of  the  RST,  each  of  the 
partial  waves  (or  normal  modes)  making  up  the  total  scattered  field 
by  a  tarqet  can  be  subdivided  into  two  contributions.  The  first  one, 
denoted  the  "background,"  represents  the  wave-field  returned  by  an 
otherwise  identical  but  impenetrable  scatterer.  The  second  portion 
of  the  scattered  wave-field  represents,  and  is  caused  by,  the 
resonances  of  the  internal  penetrable  material  making  up  the 
scatterer  (8).  The  splitting  described  above  also  holds  for  the 
reflection  coefficient  from  plane  layers  where  the  "backgrounds"  are 
usually  constants.  The  RST  was  recently  applied  to  the  problem  of 
sound  transmi ssion/reflection  through/by  a  liquid  layer  embedded 
within  one  fluid  (9),  and  later  when  the  liquid  layer  separates  two 
dissimilar  fluids  (10).  In  electromagnetism  the  method  was  first 
conceived  and  developed  for  the  spherical  geometry  (11),  and 
presently  for  sets  of  plane,  dielectric,  interacting  layers. 


I.  THE  CLASSICAL  MODEL 


We  are  particularly  concerned  here  with  the  detection  and 
possible  material  characterization  of  thin  petroleum  (oil)  layers 
lying  beneath  a  relatively  thick  overlayer  of  ice  usinq  radar  echoes. 
We  will  model  here,  however,  a  qeneral  situation  consisting  of  two 
interacting  layers  (ice  and  oil)  separating  an  upper  (air)  from  a 
lower  (water)  medium  while  (radar)  electromagnetic  waves  are  incident 
on  the  bilamlnar  system  from  the  upper  (air)  side.  (See  Figure  1.)  We 
exploit  here  the  fact  that  both  substances  have  different  dielectric 
constants  and  different  thicknesses,  which  will  excite  certain 
characteristic  "resonance  features"  in  the  returned  radar  echoes. 
These  features  can  be  made  manifest  and  much  more  understandabl e 
under  the  light  of  the  RST. 

We  assume  that  the  bilaminar  system  is  composed  of  two  ideally 
bonded  layers  with  the  upper  one  (ice)  thicker  than  the  lower  one 
(oil).  Following  the  approach  of  Brekhovskikh  (6)  the  oil-water 
interface  is  chosen  as  the  origin  of  the  coordinate  system,  as  shown 
In  Figure  1,  while  the  water  (medium  1)  and  the  air  (medium  4)  are 
considered  to  be  semi-infinite  in  the+z  directions.  Media  2  and  3, 
of  thicknesses  d2  and  d3  ,  are  the  oil  and  the  ice,  respectively,  and 
we  further  assume  that  d2 <d3 .  The  incident  electromagnetic  wave  is 
assumed  to  emanate  from  a  distant  source  In  medium  4  and  It  is 
polarized  with  the  electric  vector  £  parallel  to  the  media 


interfaces.  Maxwell  s  equations  in  Gaussian  units  are 


7  x  H 


Attq  £  +  _e_  3  E 

Co  C*  3t 


7  •  E  *  0  (1) 


7  x  H  =  - 


c0  3 1 


,  7  •  H  *  0 


where  e  is  the  dielectric  constant,  c0  is  the  speed  of  light  in 
vaccuum,  and  y,  a  are  the  permeability,  conductivity  of  the  medium 
supporting  the  wave,  respectively.  For  plane  monochromatic  waves 
with  time-dependence  e'iwt,  Maxwell's  equations  reduce  to: 


7  x  h  = 


7  x  E 


where 


e*  =  e  +  i 


In  general,  air,  fresh  water,  oil,  and  fresh  water  ice  are 
sufficiently  poor  conductors  so  that  o«.D.  We  note  that  salt  water 
and  salt  water  ice  are  relatively  better  conductors,  and  thus, 
corrections  are  needed  in  such  instances.  The  (vector)  Helmholtz 
equation  for  £  is  found  by  taking  the  curl  of  the  second  of  Eqs^. 
(2),  and  usinq  the  first  of  Eqs.  (2)  to  eliminate  the  curl  of  H. 
The  result  is 


(72  +  k2)  E 


where 


u)  /e —  ~ 

co 


Here,  w=2nf  is  the  angular  frequency  of  the  incident  wave,  and  c0  the 
speed  of  light  in  vacuum  (i.e.,  c0  =3xl010cm/sec) . 

Since  there  is  a  physical  di sconti nui ty  at  each  media  interface, 
several  boundary  conditions  must  be  satisfied.  It  is  thus  required 
that  between  the  /jth  and  the  j  +  1  media. 
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if; 


(6) 


(E.)tan  =  (E.+l)tan 

J  J 

(H.)tan  =  (H.j  +  l)tan 

J  J 


This  is  a  statement  of  continuity  of  the  tanqential  components  for  the 
electric  and  magnetic  fields  across  the  various  interfaces. 

The  electric  field  in  medium  (4)  is  given  by 


Ey  -  |a4  exp  £-1  k4(z  -  z3)cos  04J  + 
+  B4  Jexp  ik4(z  -  z3)cos  04^J|exp  £ 


iku  z  sin  6-  J 


] 


(7) 


where  and  B k  are  the  amplitude  coefficients  of  the  incidence  and 
reflected  waves,  respectively,  z,=d2  +d3  is  the  location  of  the 
interface  separating  media  (3)  and  (4),  while  04  Is  the  incidence 
angle.  The  tangential  component  of  the  magnetic  field  in  medium  (4) 
is  given,  using  Eg.  (2),  by 


H4jt  -  |  A4  exp  £-ik4(z-z3)cos  04J-  B,,exp  £ik4  (z-zj)  cos  04  J  } 

•  exp£ik4z  sin  04J 


(8) 


where 


COS04 


1 

/rr  cos04 


(9) 


is  the  field  impedance  defined  as  the  ratio  of  the  tangential 
component  of  the  electric  field  to  the  tangential  component  of  the 
magnetic  field.  Similar  expressions  are  obtained  for  the  electric 
and  magnetic  fields  in  the  other  three  media.  In  medium  (1)  one  has 
f^-0  because  there  is  no  reflected  wave.  The  connecting  bond  between 
phase-relationships  on  either  side  of  a  boundary  is  provided  by 
Snell's  law,  which  is. 
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..  t 

j 

j  knsinQu  *  kisin6i  for  i  -  1,  2,  3  (10) 


where 


fa) 

c« 


^r 


c 


Three  sets  of  equations  involving  the  coefficients  Aj  and  Bj  can 
be  obtained  usinq  the  boundary  conditions,  Eqs.  (6),  at  the 
interfaces,  [i.  e. ,  z=0,  z  =  d2,  and  z=d 2  +  d3]  In  particular,  on  the 
boundary  separating  the  2nd  and  3rd  medium,  the  boundary  conditions 
to  he  satisfied  (Eqs.  (6))  are, 


A2exp  (-Ik2d2cos02  ) 
A2exp  (-ik2d2cos  02  ) 


+  B2exp  (ik2d2eos62  ) 
-  B2exp  (lk2d2cos02  ) 


As  +  Bj 


z? 

l) 


(Aj 


-  3j) 


(12) 


Inspection  of  Figure  1  shows  that  the  reflection  coefficient  R 
is  given  by  the  ratio 


R  - 


Ah 


(13) 


Similar  expressions  can  be  obtained  for  the  continuity  conditions  at 
the  two  other  Interfaces  which  then  can  be  used  to  determine  A*  and 
Ri,  In  terms  of  Ai.  After  all  the  reductions,  the  reflection 
coefficient  of  the  bilaminar  system  under  study  is 


R  - 
where 


A~  cos  52co8  5j-B~  sin62sln6i  -  i 

rc-sindacosd,  +  D”  cos6*sin6^1 

A+  cos  52cos6|-B'r  sin  62 sin 6j  -  i 

^C+sin62co8  6|  +  D'r  cos52Sin£)J 

(H) 
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A-  «=  z2Z)(zi  +  zO 

«  zizj+  zlzi* 

^  2 
c—  *  29  (Z2  +  Z1Z4) 

D~  *  Z2  (zi  +  ZiZ),) 


and 


(15) 


62  *  “  (fdj)  £  /e2  Vl  -  ~  siJeT 

63  *  —  (fds)  /cT\fl-~^- sin20„ 

CO  £3 


(16) 


and  a  is  the  thickness  rat i 0  a  =d  3/d 2> 1 .  In  the  case  where  E  is 
parallel  to  the  air-ice  interface,  the  impedances  are 


z 


j 


cosej 


(j  -  1,  2,  3,  A) 


and  also 


(17) 


Zj 


cosdj/ /Cj 


(j  -  1,  2,  3,  A) 


for  the  case  where  E  is  parallel  to  the  plane  of  incidance. 

The  numerical  values  used  in  this  analysis  for  the  dielectric 
constants  are  as  follows: 
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Ice 


c 


Oil:  e2  =  2  ;  Water:  ei  -  81 


(18) 


In  the  linitinq  case  where  ri2-+0  (and  62  +  O),  the  reflection 
coefficient  reduces  to  that  of  a  single  layer  of  thickness  d 
separating  two  different  media,  viz, 


za(zi  -  Zu)  cosdg  -  1  tzj  -  Z1Z4.  sin  63 
z 3 (z  1  +  zi,)  cos  63  -  i  Z3  +  ZjZw  sin6j 


(19) 


This  same  result  is  also  obtained  from  Eq.  (14)  when  media  (?)  and 
(3)  are  the  same  [i.  e.  ,  z2=Z!,E!  =  c3].  For  completeness,  we  also 
give  the  transmission  coefficient  T  for  the  double  layer  system 

.  _ 2Z1Z2Z3/C0S6  2C0S6  3  _  Al 

T  “  A+  -  B+  tan  62  tan  63-  i  (C+tan62+  D+  tan  6?]  Ai,  (20) 

which  for  one  layer  (i.  e.  ,  d2-*0  62-*0)  reduces  to 


T  "  z 3 (z  1  +  zi,)  cos53-  i[z3  +  z^i^sinfis  '  (21) 

Equations  (14)  [or  (20)1  and  (19)  [or  (21)]  are  exact  solutions  for 
the  reflection  coefficient  [or  the  transmission  coefficient]  for  a 
two-layer  and  one-layer  dielectric  slab,  respectively,  separatinq  two 
semi -infinite  non-conduct i ve  media. 

The  square  of  the  magnitude  of  R  can  he  expressed  as 


1*1 


2 


zl(zi  -  zi,)*cos*<$3  +  (zi  -  zizi,)*  sin*  6j 
zj(zi  +  zm ) ^cos*  63  +  (zi  +  z  1  zi, )  *  s in' 63 


(22) 


for  a  single  layer,  and 


iRl2  - 


A~cos 62 cos 63  ~  B~sin<S;  sin  63 
’A^cosSjcos  53  -  B+sin528i-n  63J 


+  fc~sin5?cos  63  +  D~cos5?sin&;3  _ 
+  rc+sin&cos  63  -  D^cos  6zsln 63  J 
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for  the  bllaminar  conf i quration. 

Plots  of  target  cross-secti ons  or  of  reflection  coefficients 
|RJ 2  as  functions  of  frequency  (or  frequency-thickness  products)  and 
of  another  variable  such  as  angit  of  incidence  64  have  been  referred 
to  as  "response  surfaces"  in  the  past  (1?),  The  response  surface  for 
a  sinqle  (d2+0)  layer  (parallel  t  polarization)  is  shown  in  Fiqure  2 
for  frequency-thickness  products  ranqinq  from  0  to  900  MHz-meters, 
and  incidence  angles  04  from  0P  (normal)  to  90°  (qrazinq).  This 
response  surface  exhibits  a  series  of  peaks  and  valleys  which  are 
seen  to  be  evenly  spaced  in  planes  of  constant  incidence  anqles,  but 
unevenly  spaced  in  planes  of  constant  fd3  products.  The  minimum 
values  of  |R|2  correspond  to  values  for  which  the  real  part  of  the 
denominator  of  R  vanishes.  They  approximate  the  resonance 
frequencies  of  the  layer  of  thickness  d3.  For  the  case  of  a  sinqle 
layer  system,  these  resonances  occur  for  phase  delays  6s  which  are 
odd  multiples  of  tt/2.  For  normal  incidence  cases,  this  corresponds 
to  plate  thickness  which  are  proportional  to  an  odd  number  of 
quarter-wa vel enqths .  These  ei qenf requenci es  are  the  natural 
resonances  for  the  unloaded  plate. 

For  the  bilaminar  system,  the  resonances  occur  at  values  of  the 
frequency  thickness  products  fd3  and  of  the  incidence  angle  e.,  for 
which 


A+cos62Cos63  ■  B+sin52sin63  .  (24) 


For  fixed  incidence  anqles,  this  characteristic  equation 
determines  the  frequency-thickness  resonances  (fd3)3  of  the  bilaminar 
system  (1=1,  2,«..n).  Note  that  Eq.  (24)  represents  the  vanishing 
of  the  real  part  of  the  denominator  of  Eq.  (14). 

Figs.  3  and  4  exhibit  the  response  surface  for  the  bilaminar 
system,  for  various  polarizations,  ice-to-oil  thickness  ratios,  and 
in  various  frequency  regimes.  The  six  parts  of  Fig.  3  are  all 
produced  for  incidence  angles  04  in  the  range  0<94<900  and  for  the 
frequency-thickness  products  in  the  ranqe  0<fd3<  900MHz-meter .  The 
upper-most  plot  in  Fiq.  3  is  for  perpendicular  polarization,  (see 
Fiq.  1)  a  case  of  little  interest  in  this  analysis,  while  the  ot^her 
five  plots  are  all  for  the  case  of  parallel  (i.  e.,  E. .  ) 
polarization.  The  second  plot  of  Fig.  3  corresponds  to  the  case  '  'of 
a  sinqle  layer  with  parallel  polarization.  The  following  plots 
(third  to  sixth  from  the  top)  are  for  decreasing  values  of  the 
parameter  a  as  follows:  third,  for  a  =  500;  fourth,  a=100;  fifth, 
a=10;  and  sixth,  a=5.  These  correspond  to  increasingly  thicker 
layers  of  oil  (d2  )  relative  to  a  fixed  layer  of  ice  (d3).  The  last 
one,  for  example,  corresponds  to  an  ice  layer  five  tines  thicker  than 
the  oil  layer.  Figure  4  presents  again  the  response  surface,  all  for 
parallel  polarization,  but  In  the  much  larger  range  of 
frequency-thickness  products  going  from  0  to  3  G1 gaHertz-meter. 
Highly  Inclined  ridges  appear  in  the  response  surface  in  the  higher 
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frequency  parts  of  the  ranges.  The  top  plot  Is  for  a  single  layer, 
or  a  double  layer  with  a  =®.  The  second  to  fourth  plots  correspond 
to  a  values  of  a =50,000,  500  and  5,  respectively.  These  plots  are 

all  exact,  no  approximation  having  been  Introduced  yet.  Practical 

radars  operate  in  the  upper  portions  of  these  plots. 

From  Eq.  (24)  [or  Eq.  (19)],  It  follows  that  for  a  single 

layer  (i.  e. ,  62^0)  the  characteristic  equation  giving  the  resonances 

reduces  to  cosS3  =  0.  Table  presents  a  listinq  of  the  first  twenty 
roots  of  the  hilaninar  situation  governed  by  Eq.  (24).  This  table 
is  constructed  for  the  same  thickness  ratios  considered  In  Fig.  4 
which  are  a  =  d3/d2  =  50,000  ,  500,  and  5  MHz-neter.  Table  J_  is 
constructed  for  the  case  of  normal  incidence  04=0°.  Table  2  contains 
exactly  the  same  Information  as  in  Table  1  but  for  an  i nci dence  anqle 
of  04=30°.  These  entries  correspond  to  the  location  of  the  minima  in 
Fig.  4,  but  in  the  planes  04=0°  and  30°,  respectively. 


< 


II.  MODEL  FORMULATION  IN  TERMS  OF  THE  RESONANCE  SCATTER  I  NO  THEORY 

The  RST  formulation  (8)  first  linearizes  the  terms  appearinn  in 
the  rlqht  hand  side  of  Fq.  (14)  by  means  of  a  Taylor  series 
truncated  after  Its  second  term.  We  denote 

Y~  (fdi)  “  A"  cosSicosfis”  B*  sinSjsinfij  (25) 


then,  the  two-term  Taylor  series  expansion  about  fQn  is  (for  each  n): 


F*  (fd,)  -  F* 


3P+  1 

(fond»>  +  (f"  fon>  ds§(fd7)Jt  r“ 

*“*on 


(26) 


Here,  fon  are  The  resonance  frequencies,  roots  of  Eq.  (24)  for 
constant  incidence  anqle  04  ,  or  which  are  the  same,  they  are  the 
roots  of 


F+(fds)  -  0  ,  (for  f  -  f0n>  <27> 


It  is  now  clear  that 


F"(fonds)  *  A“cos  62o'^cos63?^  -  B”sin6  »^sin5^so\  (28) 


and 


3F-(fonds) 

2(f0nd*) 


‘  -  “*,+ 


V4  3(fds)  +  ^fdal^ 


sin52o^  cosfiso^ 


] 


Lon 


(29) 


and  the  derivatives  are  found  from  Eqs.  (16)  to  he: 


4fH - ?!  ±  /ZT  <1  -  .ln!0»  )■» 

B(fds)  co  a  C2 


(30) 


and 


«Kfdj)  co 


1 3 


(31) 
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It  is  not  hard  to  show  that  Eq.  (14)  can  now  be  re-written  as 


where 


^  a  +  b(x  ~  xon)  "*■  *  s(^n/2) 

Rm  2^/  x  -  xon  +  i  (rn/2) 

n*l 

~f  m  ~  f  sV+l  [c+sin6a>n)cos6jjn)  +  D+cos6a»n)sin6^n^J 
.  3  x  .  xon 


(32) 


(33) 


b  ** 


(34) 


a  C~  sln(.21p^cos6^  +  D*‘cos  62^sini5sf^ 
C4  sin^cosfil?^  +  D+cos  S^sinfi^ 


and  x=fd3. 


(A~B+  -  A^B-jsinS^^inS^ 


(36) 


The  form  of  the  reflection  coefficient  in  Eq.  (32)  exhibits  the 
resonance  character  of  the  phenomenon,  it  was  borrowed  from  nuclear 
scattering  theory,  where  similar  expressions  a  e  associated  with  the 
names  of  Rreit  and  Wiqner  (13),  and  it  is  now  the  central-type  of 
expression  for  the  RST-analysis  of  macroscopic  scattering  problems  in 
electromagnetics,  acoustics,  and  el  a stodynamics  (8).  The  summation 
in  Eq.  (32)  is  merely  symbolic,  since  the  formula  and  the 
RST-approximati on  are  only  valid  in  the  vicinity  of  the  resonances, 
and  the  sum  merely  means  that  the  neighborhoods  of  all  the  resonances 
must  be  aciounted  for.  Each  resonance  xon  has  a  width  rn  associated 
with  the  imaginary  part  of  the  simple  denominator  of  Eo.  (32). 


In  terms  of  the  RST  formalism,  it  follows  from  the  above 
analysis,  that,  the  square  of  the  modulus  of  R  which  is  the  quantity 
plotted  in  all  the  graphs,  is  given  by 
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a  ,+.  Mx  -  x^n)  I 

(x  -  *.nn)  +  rT 


for  the  bllamlnar  system.  For  the  single  layer,  this  expression 
reduces  to 


e2 (x  -  x0n)2  +  f*  Yr 
<x  "  xon>  +  Yn */4 


where  the  x  on  are  now  the  roots  of  cos63=0,  rather  than  the  roots  of 
the  characteristic  equation  (?4),  and  where  a-*T)  and 


e  «  (zi  -  z4)/(zi  +  z4 )  =  A“/A+ 

f  =  (zl  -  z;z4)/(zl  +  2jz4)  =  D~/D"* 


Z3(Z1  +  Zu) 


sin204 
^  3 


and  the  layer  of  thickness  d2  (i.  e.  ,  oil)  has  been  removed  since  in 
this  case  d2-*0.  We  have  already  shown  in  the  acoustic  case  (14)  how 
the  casting  of  the  expressions  into  the  RST-forms  can  he  made  exact 
hy  means  of  Mi ttao-Lef f 1 er  rather  than  Taylor  series  expansions,  a 
point  that  will  not  be  considered  further  here,  since  the  approximate 
RST  expressions  developed  here  are  adequately  qood  and  useful  for  our 
purposes.  We  note  that  the  one-layer  resonance  relation  in  Eq. 
(38),  obtained  as  a  particular  case  of  our  more  general  Eq.  (32),  is 
the  electromagnetic  counterpart  of  an  available  acoustic 
relation  (1°)  useful  to  study  ocean  bottom  properties. 

III.  NUMERICAL  CALCULATIONS  AND  DISCUSSION 

The  quantity  JR)2  is  a  scalar  function  of  incidence  anqle  84  and 
frequency-thickness  products  fd3.  This  quantity  can  be  plotted 
three-dimensional ly  as  we  have  shown  before  in  Figs.  2,  3,  and  4. 
It  can  also  be  plotted  two-dimensional ly  versus  one  of  the  variables, 
keeping  the  other  one  fixed.  For  normal  Incidence  (i.  e.,  84  =0), 
Fig.  5  shows  the  two-dimensional  response  curve  (i.  e. ,  |R|  vs. 


\ist. 
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fd3)  in  the  range  0<fd<  1,000  megahertz-meters,  first  for  a  single 
layer  of  ice  separating  air  from  water  (and  for  Pii  polarization), 
and  later  for  a  double  layer  of  ice  and  oil  with  four  levels  of 
oi 1 -t h i ck nes ses  relative  to  the  Ice's.  The  relevant  values  of  o  In 
the  bottom  four  graphs  of  Fig.  5  are  ctsdj/d2"500,  100,  25,  5, 

respect i vely .  Two  sets  of  curves  are  shown  In  these  graphs,  one  in 
solid  and  the  other  in  dashed  lines.  The  solid  lines  represent  the 
results  of  exact  calculations  using  the  (complicated)  classical 
formulism  that  led  to  Eqs.  (23),  (15),  and  (16).  It  is  quite 

difficult  to  interpret  the  physical  phenomenon  from  these  cumbersome 
expressions.  The  dashed  lines  represent  the  results  of  approximate 
calculations  using  the  RST-formul 1 sm  that  led  to  Eqs.  (33)  -  (37). 
It  is  quite  evident  from  these  plots  that  the  RST-formul 1 sm  generates 
each  resonance  "dip",  one  at  a  time,  and  that  the  agreement  between 
dashed  and  solid  curves  near  the  tip  of  each  minima  (or  dip)  Is 
nearly  perfect,  and  only  away  from  each  resonance  minima  do 
differences  become  apparent  between  the  two.  In  fact,  the  agreement 
betwen  solid  and  dashed  curves  is  near  perfect  down  to  the  half-width 
of  each  resonance-dip,  measured  at  half-minimum,  and  given  by  Eg. 
(33). 

The  RST-form  in  Eq.  (37)  or  better  in  (32),  is  much  clearer  to 
understand  and  it  represents  the  addition  of  cont ri but i ons  each 
emerginq  from  a  pole  located  at 

x  3  xon  “  1  rn/2  (40) 

in  the  complex  x-plane.  We  have  expanded  on  this  type  of 
interpretation  before  (15)  in  similar  Instances.  The  real  part  of 
the  pole-position  qives  the  location  of  the  resonance  minimum  In  the 
x=fd3  axis,  and  its  imaginary  part  is  its  half-width.  As  shown  in 
Fia.  5  as  well  as  in  Table  1,  for  a  single  layer  the  spacing  between 
resonances,  the  magnitudes  of  |R|2  at  resonance,  and  the  resonance 
half-widths  are  all  constants.  Further  observation  of  Fig.  5  shows 
that  when  the  ice  layer  is  very  thick  compared  to  the  oil  layer 
(1.  e.,a=500),  there  is  hardly  any  visible  difference  between  the 
response  curve  for  the  bilaminar  configuration  and  that  of  the  single 
ice  layer  shown  at  the  top  of  Fig.  5.  The  resonance  locations, 
roots  of  Eq.  (24),  approach  those  of  the  single  dielectric  slab  as 
d2-*  0.  For  ice-to-oll  thickness  ratios  of  100  or  less,  the  response 
curve  of  the  bilaminar  layer  configuration  starts  to  differ  visibly 
from  that  of  the  single  layer.  The  difference  manifests  Itself  in  an 
upward  slop*  of  the  response  curve  that  culminates  In  the  wavy 
undulations  shown  at  the  bottom  of  Fig.  5,  where  a*5.  The  pattern 
exhibited  is  similar  to  the  modulation  produced  by  the  beating 
toqether  of  two  sinusoidal  waves  of  different  frequencies.  We  repeat 
that  all  results  displayed  In  Fig.  5  exhibit  excellent  agreement 
between  the  exact  solution  (solid  lines)  and  the  approximations 
inherent  to  the  Resonance  Scattering  Theory  (dashed  lines).  We 
remark  that  there  is  a  minute  difference  between  the  exact 
resonances,  at  the  minima  points  of  the  solid  line  graphs,  and  the 
RST  (dashed)  approximations,  which  Is  due  to  the  presence  of  the 
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quantity  "a"  appearing  In  the  numerator  of  Eq.  (32),  which  is  small 
for  the  case  at  hand.  We  also  note  that  while  the  spacing  between 
resonances  and  their  widths  were  constants  for  the  single-layer  case, 
this  no  longer  holds  for  the  bllaminar  configuration  where  a  "cyclic" 
variation  is  now  present. 


IV.  THE  INVERSE  SCATTERING  PROBLEM 


It  is  assumed  that  the  electromagnetic  properties  of  air  and 
water,  in  this  case  £4  ,  £i  ,  are  known.  The  solution  of  the  inverse 
scattering  problem  then  reduces  to  a  determination  of  the  paramters 
d2 ,  £2  ,  d  3,  £3  for  the  thicknesses  and  dielectric  constants  of  the 
oil  and  ice  layers.  We  will  show  that  sufficient  information  can  be 
extracted  from  a  careful  analysis  of  the  response  surfaces  and  curves 
presented  in  Figs.  2  through  6  to  completely  identify  these 
parameters.  The  key  to  the  solutionlies  in  an  understanding  of  the 
characteristic  Eq.  (24)  which  determines  the  locations  where 
resonances  occur.  The  quantities  A+  and  B+  present  in  this  equation 
depend  on  the  dielectric  constants  of  the  four  media  involved  in  the 
reflection  process,  as  well  as  on  the  initial  incidence  angle  e4 . 
For  simplicity  in  the  anlaysis  we  shall  assume  normal  incidence 
although  any  oblique  incidence  could  have  been  used  just  as  well. 
The  phase  delays  introduced  by  the  two  finite  layers,  62  and  63,  (c/o 
Eqs.  16)  depend  on  the  dielectric  constant  and  thickness  of  the 
respective  layer  as  well  as  on  the  frequency  of  the  incident 
radiation.  We  note  that  if  the  thickness  of  the  oil  layer  is  small 
and  the  frequency  is  low,  then  a,  will  have  a  small  value.  If  it  is 
further  required  that  63  not  be  small,  then  the  characteristic  Eq. 
(24)  is  approximated  by 

cosd3  *  0  (41) 


which  we  saw  was  the  resonance  condition  or  characteristic  equation 
for  a  single  (upper)  layer  (of  ice).  Thus,  at  low  frequencies,  the 
hilaminar  layer  problem  with  a  large  ratio  a  of  ice-to-oil 
thicknesses,  behaves  to  first  order  like  a  single  ice  layer. 
However,  as  the  frequency  increases  or  as  the  oil  layer  thickness 
increases,  the  effects  of  the  oil  layer  on  the  resulting  response 
curves  become  more  pronounced.  There  is  a  modulation  effect  visible 
in  the  response  curves  (c/o  Fig.  5)  and  caused  by  the  presence  of 
the  oil  layer  which  was  discussed  in  §1 1 1 .  We  further  note  that  in 
the  bottom  graph  of  Fig.  5,  the  first  displayed  resonance  is  very 
near  in  amplitude  to  that  of  the  single  ice  layer  given  in  the  top 
graph  of  Fig.  5.  As  the  frequency  increases,  the  amplitude  of  the 
following  resonances  start  to  increase  approaching  a  minimum  which 
coincides  with  the  value  seen  for  a  single  oil  layer  in  Fig.  6.  How 
closely  the  resonances  will  approach  those  of  the  single  oil  layer  is 
ultimately  governed  by  Eq.  (24).  It,  therefore,  followsthat  the 
amplitude  oscillations  introduced  by  the  modulation  effect  here 
described  (c/o  Fig.  5,  bottom  graph)  are  bounded  by  the  extreme 
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values  of  the  resonance  amplitudes  of  the  oil  or  the  ice  layers 
individually  and  the  modulation  Is  confined  to  occur  in  between  them. 

It  can  be  seen  from  Eq.  (32)  that  at  any  resonance 


R 


4rF' 


[(rn/2)a'+  8*] 


(42) 


where  1^/2,  q,  and  a  were  given  be  Eqs.  (33),  (35)  and  (36), 
respectively.  It  is  easily  shown  for  resonances  where  either  cos  62 
or  cos  63  is  close  to  zero,  that  a/(Fn/2)  is  negligibly  small.  By 
means  of  Eqs.  (29)  -  (36)  this  ratio  Is 

a  _ _ (A~B+  -  A+B")  cos6  2Q^cos  63^ 

(rn/2)  B+  (c+sin^sind^V  D+cosd^  ain&yo) 

A+B~  -  A~B+  .  _ 1 _ 

[c+tan  D+tat6jo^l 

Thus,  for  large  values  of  either  tanfi^  or  tand^jl,  it.  follows  that 


|a/(rn/2)l  <  [tan  or  [tan  (44) 


Hence,  In  the  neighborhood  of  these  resonances,  the  ratio  |  2a/Fn  |  is 
negligible  and  any  resonance  amplitude  is  essentially  given  by 

R  a  8 


which  reduces  to 


Ro  * 


(45) 


for  an  oil-type  resonance 


provided  tha\  tan62>>l,  and  to 


R  .  _21_  (46) 

Ri  D+ 


for  an  ice-type  resonance  provided  that  tan6j>>l. 
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For  the  specific  problem  considered  here  (i.  e.  ,  ice  and  oil) 
the  resonances  with  the  smallest  amplitudes  (i.  e.,  |  R  |2  ==0.‘226)  are 
caused  by  the  ice  layer  (see  Fig.  5  top)  while,  the  resonances  with 
the  laroest  amplitudes  (i.  e.,  |R|2-  0.405)  are  due  to  the  oil  layer 
(See  Fiq.  6).  If  the  amplitudes  of  the  ice  and  oil  resonances, 
distinguished  in  this  fashion,  are  represented  by  Po  and  Ri, 
respeci  t  i  vely ,  then  it  follows  from  Eqs'.  (45)  and  tire  expressions 
for  C  +  or  D±  in  Fqs.  ( 1 5  )  t  h  a  t. 

..  -  <*.“>■»  (-Krr)15  •  (47) 

The  impedances  2\  and  2  <,  are  known  quantities,  Fq.  (9),  while 
I'o  is  determined  from  the  calculated  response  surfaces  of  from 
patterns  similar  to  those  in  Figs.  5  and  6  obtained  from 
measurements.  For  normal  incidence,  Eq.  (47)  directly  yields  the 
dielectric  constant  of  the  hidden  (oil)  layer  in  the  form 


£2  =  Uiz-.)-1* 


(48) 


In  a  similar  fashion,  the  dielectric  constant,  for  ,the  upper 
(ice)  layer  is 


(49) 


In  what  follows  we  will  assume  these  dielectric  constants  known 
and  niven  by  Eqs.  (4P),  (49).  Thus,  all  the  constants  A±,  P±.  C±, 
D+  appearing  in  Fq.  (15)  can  be  determined  from  selected  resonance 
amplitudes  in  the  reflection  coefficient  plots.  The  rema i ni no  still 
undetermined  quantities  are  the  thicknesses  of  the  oil  and  the  ice 
layers.  This  information  can  be  extracted  from  the  (half)  widths  of 
the  resonances. 

For  the  purposes  of  solvinq  the  inverse  scattering  problem,  we 
note  that  the  response  surfaces  and  curves  in  Figs.  2-6  are  all 
really  found  in  practice  as  functions  of  frequency  if,  rather  than 
frequency-thickness  products  fd3.  Figs.  2-6  have  all  been  displayed 
versus  x  ( =fd3  )  for  convenience  In  the  formulation  and  ease  in  the 
display,  but  with  the  understandi ng  that  f  was  the  physical  quantity 
varying  alonq  those  axes  since  d3  is  not  really  known  apriori.  The 
difference  between  x  and  f  is  the  constant  scale  factor  d3.  If  those 
response  curves  were  displayed  versus  f,  the  amplitudes  of  the 
resonances  would  remain  the  same  as  in  Figs.  2-6  but  their  width 
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would  be  changed  by  a  factor  of  d3.  The  widths  Tn  g.iven  iq  Eq.  (33) 
in  plots  displayed  versus  x  will  now  become  widths  in  plots 
displayed  versus  f  and  they  are  given  by 


Tri  m  _  C^sinS^  cos  +  D^cos  620^  sin 

r  a  F+i  ’ 

d  j  — « — 

3x 


(SO) 


'on 


a  value  which  differs  from  that  in  Eq.  (33)  by  a  scale  factor  d3 
present  in  the  denotni  nator .  It  is  clear  that 


ds 


3F+ 

3x 


and  that 


|[a+  /t,  c-‘  +  b+/ei]  stosSs’  cosfS;1 

+  B+  /eTa“iJ:os  620  sin53o| 


+ 


(51) 


For  an  oil-type  resonance,  for  which  tan<5i>>l,  Eq.  (50)  gives 
approximate  result 


jio 


2 


(52) 


and  for  an  ice-type  resonance,  for  which  tan63>>l, 


now 


the 


Dividing  these  two  equations  one  by  the  other  and  solving  the  result 
for  a-1  (=d2/d»)  yields, 


where 


A+S  -  B+ 

tr  -  sir 


s  - 


c* 

tF 


ft 


(54) 


(55) 


Once  a’1  is  known,  either  Eq.  (52)  or  (54)  can  be  used  to  determine 
d  2  and  d3 .  These  thicknesses  come  out  to  he 


da* 


dj- 


(56a) 


(56b) 


which  are  expressed  here  In  terns  of  measured  widths  r„#  ,  rni  and  of 
quanltltes  such  as  A  +  ,  B ♦,  C+,  D+,  gz  ,  e,  ,  all  determined  earlier. 
All  the  material  and  geometrical  parameters  required  to  completely 
solve  the  Inverse  scattering  problem  have  been  determined  from  the 
radar  reflections. 
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would  be  changed  by  a  factor  of  d3.  The  widths  qiven  in  Eq.  (33) 
in  plots  displayed  versus  x  will  now  become  wi dths  r ^  in  plots 
displayed  versus  f  and  they  are  qiven  by 


a  value  which  differs  from  that  in  Eq.  (33)  by  a  scale  factor  d3  now 
present  in  the  denominator.  It  is  clear  that 


and  that 


3F+  \  _  _2tt 

3  x  7  co 

k  '  X _ 


-  {[a+  VT2  a-1  +  B+  /El]  Sin 5^  cosfft* 

n  '  . 

£a+  /ti  +  B+  /TTct-Jos  620  sin  5jo| 


For  an  oil-type  resonance,  for  which  tan<$j>>l,  Eq.  (50)  gives  the 
approximate  result 


■Jkkm 

2 


and  for  an  ice-type  resonance,  for  which  tan6a>>l. 


_ _ 

-f-  d,  [a+/ET  +  b+  StT  cr1] 
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Dividing  these  two  equations  one  by  the  other  and  solving  the  result 
for  or1  (=d2/dj)  yields, 


where 


(54) 


(55) 


Once  a"1  is  known,  either  Eq.  (52)  or  (54)  can  be  used  to  determine 
d  2  and  d3 .  These  thicknesses  cone  out  to  be 


d3* 


(56a) 


(56b) 


which  are  expressed  here  In  terns  of  measured  widths  Tno  ,  rni  and  of 
quanltites  such  as  A*,  B +,  C+,  D\  e2  ,  e,  ,  all  determined  earlier. 
All  the  material  and  geometrical  parameters  required  to  completely 
solve  the  inverse  scattering  problem  have  been  determined  from  the 
radar  reflections. 


V.  CONCLUSIONS 


We  have  constructed  the  RST  predictions  for  the  reflection  coefficient 
of  a  bllamlnar  plane  configuration  of  dielectric  layers  separating  two 
dissimilar,  non-conducting  media.  This  predlcltlon  Is  compared  to  the 
classical  model  prediction.  The  comparlslon  serves  to  emphasize  the  accu¬ 
racy  and  simplicity  of  the  Resonance  Scattering  Theory  (8),  as  well  as  Its 
ability  to  yelld  a  clear  physical  Interpretation  of  a  complex  phenomenmom  In 
terms  of  scattering  poles.  Furthermore,  the  RST  shows  Its  capability  as  a 
material  probe  to  extract  geometrical  and  material  composition  Information 
from  the  set  of  scattering  layers  returning  the  reflected  radar  echoes. 
This  capability  is  used  to  completely  solve  the  Inverse  scattering  problem 
for  the  bllamlnar  configuration  under  analysis.  The  locations  and  widths  of 
the  resonances  as  well  as  their  amplitudes  determine  the  thickness  and 
dielectric  compositions  of  both  Interacting  layers.  The  process  disen¬ 
tangles  which  set  of  resonance  features  in  the  reflected  echo  Is  caused  by 
which  of  the  layers.  The  calculations  that  we  display  In  many  graphs 
Illustrate  the  case  of  an  oil  layer  hidden  under  a  covering  of  ice,  both 
media  sandwiched  between  air  on  top  and  water  at  the  bottom.  This  work 
could  be  useful  In  the  further  development  of  techniques  for  detecting  and 
mapping  the  extent  of  an  oil  spill  under  solid  Ice. 
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TABLE  1 


CASE;  Normal  Incidence  <i.e.,  64  >0°) 

LOCATION  OF  THE  RESONANCE  FREQUENCIES 


SINGLE  LAYER  DOUBLE  LAYER 


l 

(MHz-meters) 

a  »  50,000 

a  -  500 

a  «  5 

I 

41.9 

41.9 

41.9 

35.2 

2 

125.8 

125.8 

125.6 

106.4 

3 

209.6 

209.6 

209.3 

179.2 

4 

293.5 

293.5 

292.9 

253.1 

5 

377.3 

377.3 

376.6 

372.2 

6 

461.2 

461.2 

460.3 

400.5 

7 

545.0 

545.0 

544.0 

472.2 

6 

628.9 

628.9 

627.9 

542.7 

9 

712.8 

712.7 

711.4 

613.5 

10 

796.6 

796.6 

795.1 

685.7 

11 

880.5 

880.4 

878.8 

759.4 

12 

964.3 

964.3 

962.5 

833.6 

13 

1048.2 

1048.1 

1046.1 

907.2 

14 

1132.0 

1132.0 

1129.8 

979.4 

IS 

1215.9 

1215.8 

1213.5 

1050.2 

16 

1299.7 

1299.7 

1297.2 

1120.7 

17 

1383.6 

1383.5 

1380.9 

1192.4 

18 

1467.4 

1467.4 

1464.6 

1265.7 

19 

1551.3 

1551.2 

1548.3 

1339.8 

20 

1635.1 

1635.1 

1631.9 

1413.7 

22 


TABLE  2 


CASE:  Oblique  Incidence  at  6i*  *  30° 


LOCATION  OF  THE  RESONANCE  FREQUENCIES 


SINGLE  LAYER 

DOUBLE  LAYER 

£ 

(MHz-meters) 

a  -  50,000 

a  =  500 

a  «  5 

1 

43.7 

43.7 

43.6 

36.7 

2 

131.0 

131.0 

130.7 

110.9 

3 

218.3 

218.3 

217.9 

187.0 

4 

305.7 

305.7 

305.1 

264.4 

5 

393.0 

393.0 

392.2 

342.0 

6 

480.3 

480.3 

479.4 

418.9 

7 

567.7 

567.7 

566.6 

494.0 

8 

655.0 

655.0 

653.7 

567.5 

9 

742.3 

742.3 

740.9 

641.1 

10 

829.7 

829.7 

828.1 

716.3 

11 

917.0 

917.0 

915.2 

793.1 

12 

1004.3 

1004.3 

1002.4 

870.8 

13 

1091.7 

1091.6 

1089.6 

948.2 

14 

1179.0 

1179.0 

1176.7 

1024.2 

15 

1266.3 

1266.3 

1263.9 

1098.5 

16 

1353.7 

1353.6 

1351.1 

1171.8 

17 

1441.0 

1441.0 

1438.2 

1246VO 

18 

1528.3 

1528.3 

1525.4 

1322.1 

19 

1615.7 

1615.6 

1612.6 

1399.5 

20 

1703.0 

1703.0 

1699.8 

1477.2 

23 


SINGLE  AND  BI LAMINAR  CONE I GURAT IONS  FOR  FREQUENCY  THICKNESS  PRODUCTS  TO  3000  MHZ -METERS 
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FREQUENCY  THICKNESS  (MHz-M) 

RESPONSE  CURVES  FOR  SINGLE  AND  B I LAMINAR  CONFIGURATIONS 
SOLID  LINE  FOR  CLASSICAL  THEORY,  DASHED  LINE  FOR  RST. 
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FREQUENCY  THICKNESS  (MHZ-M) 

RESPONSE  CURVE  FOR  SINGLE  LAYER.  SOLID  LINE  FOR  CLASSICAL  THEORY,  DASHED  LINE  FOR  RST. 


